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We study inelastic (dynamical) impurity scattering effects in two-band superconductors with 
the same (s++ wave) or different (s± wave) sign order parameters. We focus on the enhancement 
of the superconducting transition temperature T c by magnetic interband scattering with the 
interchange of crystal-field singlet ground and multiplet excited states. Either the s++-wave 
or s±-wave state is favored by the impurity-mediated pairing, which depends on the magnetic 
and nonmagnetic scattering strengths derived from the hybridization of the impurity states 
with the conduction bands. The details are examined for the singlet-triplet configuration that 
is suggestive of Pr impurities in the skutterudite superconductor LaOs4Sbi2. 
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1. Introduction 

The study of impurity effects on multiband super- 
conductivity is stimulated by the discovery of high- 
temperature (high T c ) superconductors with FeAs lay- 
ers, 1 ~ 8 ) while it has been performed for another high- 
T c superconductor, MgB2. 9 ~ n ) One of the interesting 
points for the impurity problem is their anisotropic band 
property. In particular, much attention is paid to the 
s_|--wave pairing state characterized by sign reversal of 
order parameters with full gaps. 12-16 ) For a depairing 
effect caused by impurity scattering, the s±-wave state 
behaves like d-wave pairing rather than the conventional 
s-wave state. 17 ) Multiband properties are also reported 
in the heavy-fermion superconductor PrOs4Sbi2 whose 
higher T c than LaOs4Sbi2 implies a crucial role of the 
Pr /-electron states. 18 ' 19 ) Interband scattering due to 
impurities may be important in the multiband, the roles 
of which remain to be elucidated. 

Recently, Senga and Kontani have studied the effects 
of intraband and interband nonmagnetic impurity scat- 
terings in the s^-wave state using a simple two-band 
BCS model. 20 ' 21) They found that T c is not markedly 
suppressed when the intraband and interband scatter- 
ing strengths are not equal and are sufficiently strong. 
We also applied a similar model in a magnetic impu- 
rity case and solved a single-impurity problem. 22 ) For 
the interband scattering in the s±-wave state, the roles 
of magnetic impurities are equivalent to those of nonmag- 
netic impurities for the intraband scattering. 17 ) Very re- 
cently, it has been reported that T c is not markedly sup- 
pressed so much by the interband magnetic scattering. 23 ) 
These results imply a possibility of T c enhancement 
due to magnetic impurities having an internal structure. 
Such enhancement was argued by Fulde et al. for non- 
magnetic impurities in single-band superconductors. 24 ) 
An attractive interaction stems from inelastic nonmag- 
netic scattering by impurities, which is analogous to the 
electron-phonon origin of BCS superconductors. T c can 



be enhanced by doping such impurities with appropriate 
crystal-field level splitting, while it is suppressed by mag- 
netic impurities. 24 ) This idea has recently been applied to 
a skutterudite superconductor Lai_ a; Pr a ;Os4Sbi2 to ac- 
count for the T c enhancement by Pr substitution for La 
in LaOs4Sbi2- 25 ) Since only a single band is considered 
there, one always finds that magnetic impurities cause 
T c suppression. However, this conventional understand- 
ing has to be reexamined for multiband superconduc- 
tors. In the case of s±-wave superconductivity, magnetic 
impurities for intraband scattering behave like nonmag- 
netic impurities for interband scattering, and vice versa. 
Accordingly, we can expect T c to be enhanced by the 
inelastic (dynamical) magnetic scattering. 

In this paper, we extend the theory of Fulde et al. 24 ' 
and show a case of magnetic impurities in the s±-wave 
state in §2. Considering spin-dependent interband scat- 
tering, we examine a self-energy in the Born approxi- 
mation to derive a gap equation for two bands. The in- 
crease in T c is associated with the sign change of the 
superconducting order parameter via the inelastic impu- 
rity scattering. It is demonstrated for the s±-wave state 
whose T c is enhanced by the magnetic interband scatter- 
ing due to the singlet- singlet configuration. In §3, we also 
show some typical examples of the corresponding inter- 
band scattering. There are various octupolar (combina- 
tion of spin and orbital) scattering types that increase 
T c accompanied by an interchange between two crystal- 
field singlet states. The same argument is applied to the 
singlct-multiplet configuration that can be realized for 
such /-electron impurities as rare-earth or actinide ions 
embedded in cubic or uniaxial anisotropic (tetragonal or 
hexagonal) crystals. For the singlet-doublet, we discuss 
the s-t-wave and s++-wave combination via intraband 
scattering. The s ++ -wave pairing is characterized by the 
same sign order parameters of the two bands. We show a 
possible crystal-field splitting for T c enhancement, how- 
ever, T c is suppressed by the competition between mag- 
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netic and nonmagnetic scattering effects. For the singlet- 
triplet, we find a case in which the s ++ wave is favorable 
for the T c enhancement by magnetic intcrband scatter- 
ing. Which is chosen, the s± wave or the s++ wave, by the 
impurity- mediated pairing depends on (1) the scatter- 
ing type (dipole, quadrupole, octupole, etc.) and (2)thc 
hybridization between the impurity atomic orbitals and 
conduction bands. The details are described for an impu- 
rity with the singlet-triplet configuration. Its application 
to Pr impurity effects in the LaOs4Sb 12 superconductor 
is discussed in §4. Conclusions are given in §5. 

2. Formulation for T c Enhancement 

In this section, we present a formulation for T c en- 
hancement by inelastic impurity scattering in multiband 
superconductors. First, we briefly review a work by Fulde 
et al. for single-band superconductors. 24 -' Then, we ex- 
tend it to a two-band system as a simple case of multi- 
band and give an example to understand multiband ef- 
fects on T c . 

2.1 Inelastic impurity scattering and T c enhancement 
in single-band superconductors 
Let us begin with the model Hamiltonian % = %c + 
'Hi+'H' that consists of conduction electron "He , impurity 
"Hi, and impurity scattering W terms. The first term is 
written as 



U c = 






dripl(r)e(-iW)ip a (r) 



A / dr 



4(r)V|(r)+^(r)V t (r-)l. (1) 



Here, ipcr(r) is a field operator of the conduction election 
of the cr(=t, 4-) spin whose kinetic energy is expressed by 
the operator e(iV) = — V 2 /2m c — E-p measured from the 
Fermi energy Ep , where m c represents the electron mass 
and the Planck constant h is taken as unity. A is the s- 
wave superconducting order parameter that we assume 
to be a real value. The Hamiltonian %i for the impurity 
states is given by 



Hi 



¥ E " 



mQ"yrn®"~f'm' 



(2) 



Here, Ry represents the position of the 7th impurity. 
ajj, m and a im are the pseudo-fermion creation and an- 
nihilation operators, respectively, for the mth impurity 
energy level S m at the 7th impurity site. 26 ) The interac- 
tion Hamiltonian at the impurities is defined by 

W = J2J2Y.J dra\ m a 7n 6(r - J R 7 )M„ m ^(r)V'a(r). 



R~ ran a 



(3) 



We consider here only nonmagnetic impurity scattering 
since magnetic impurity scattering cannot enhance T c in 
a single-band case. 24 - 1 M mn is a matrix element that de- 
scribes the scattering of conduction electrons accompa- 
nied by an interchange among the mth and nth energy 
levels. 

We introduce the following 4x4 matrix form of the 



thermal Green's function: 

G{T,r,r') = -(T*(r,T)* f (r',0)), 



(4) 



where \&(r) and M*' (r) arc four-dimensional vectors de- 
fined as 



¥(r) 



/ Mr) \ 

M r ) 

4(r) 



* f (r) - ( 4(r) i>\(r) Mr) Mr) ) , (5) 
with their Heisenberg representations 



*(r,r)=e HT *(r)e 



-Ht 



^{r,T)=c HT ^{r)c- HT 
(6) 

In the absence of impurity scattering, the unperturbed 
Green's function is Fourier-transformed to 

\uji + e k p 3 + Ap 2 d- 2 



G {m,k) = -- 



A 2 



(7) 



where a a and p a (a is denoted by 1, 2, and 3 instead of x, 
y, and z, respectively, hereafter) are the Pauli matrices 
for the spin space and particle-hole space, respectively. 
Similarly, the matrix for impurity scattering is given by 



U ran — ^rt 



iP3- 



(8) 



Following Fulde et al., 24 ' we study the T c enhancement 
on the basis of the second Born approximation. For the 
scattering matrix U mn in eq. (8), Fig. 1 shows the self- 
energy given by 

t(iui) = -n imp T 2 ^ ^2 



ran uj\uj2 



IWl - #m 1^2 - $n 



k 



U mn G (iuji + iwi - iw 2 , k)U n 



(9) 



Here, 7ii mp represents the impurity density. l/(iwi — 5 m ) 
is an unperturbed Green's function for the mth impu- 
rity energy level. u>i and L02 are Matsubara frequencies 
for fermions. f2 represents the system volume. The Boltz- 
mann constant /cb is taken as unity. The important point 
is that A in eq. (7) changes its sign with the U mn GoU nm 
transformation. This is a key to the T c enhancement due 
to the inelastic impurity scattering, by analogy with the 
optical phonon case leading to an attractive interaction 
for pairing. 24 ' On the other hand, for the elastic impurity 
scattering, the pairing interaction is cancelled out by a 
depairing effect. 

In the presence of the impurities, the linearized gap 
equation is given by 



Alog-^=.T c ]TA 



S A (iw;) -i 



.^(iwi) 



to, 



(10) 



where T c (T c q) is the transition temperature in the pres- 
ence (absence) of impurities. Ea and S u are self-energies 
corresponding to the order parameter and Matsubara fre- 
quency components, respectively. 24 ' In Green's function, 
the rcnormalizcd frequency uji and the order parameter 
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CO, 




and to that for the magnetic scattering (jvi is denoted 

by r 12 ), 



— T"l2 lo S TFT" = /a(x) + Mx) 
T J-cO 



■f.ix) < 0. (15) 



co 2 -co 



CO, 



co,-co 7 + co 



'1 



CO, 



Fig. 1. Fcynman diagram of self-energy in the second Born ap- 
proximation. 



A; are given by 

lui = ui + iS 



(iwz), 



A, = A[l + £ A (iw,)]- (11) 



The self-energy in Fig. 1 leads to different signs of £a 
in eq. (10) between time-reversal invariant scattering 
and non-time-reversal invariant scattering. 24 > The former 
(nonmagnetic scattering in a single band) contributes to 
T c enhancement. In a simple case of two energy levels 
(62 > Si), after the lengthy calculation of the self-energy, 
we obtain each term on the right-hand side of eq. (10) as 

— rS A (iu;;) = -— /a (a;), 



■kT c 



E 
1 



Ul 



8T c r 12 



tanhx 



^e 



i S w (icj;) 



1 



l 



M 



CO, 



-A(x)--B(x), 



(12) 



8T C T1 2 



f u (x) = -1+tanh 2 



\b{x), 



These equations show that for any finite crystal-field 
splitting, T c is enhanced by nonmagnetic impurities in 
single-band s-wave superconductors, while it is sup- 
pressed by magnetic impurities. 24 ^ For the nonmagnetic 
scattering in eq. (14), we note that T c keeps T c0 for x = 
(elastic scattering) that gives /o(0) = [/a(0) = fu,{0) = 
-1], while \f A (x)\ > \f u (x)\ [/a (x) < 0,f w (x) < 0] for 
x > leads to T c > T c0 . Once the two levels are split, 
the exchange scattering process gives rise to an attrac- 
tive interaction for pairing, by analogy with the optical 
phonon effects in superconductivity. 24 ) 

We would like to mention a linear combination of 
cqs. (14) and (15), [fo(x)/r® 2 + / s (^)/t 12 ], which is de- 
rived from the singlet-multiplet configuration since both 
magnetic and nonmagnetic scatterings coexist. In this 
case, the weight of /a(x) becomes smaller in the com- 
bined equation, so that x has to be sufficiently large 
to satisfy the condition |/a(£)| > |/ w (a;)| for T c en- 
hancement. The details will be examined for the singlet- 
doublet configuration in §3.2. 

2.2 Extension to multiband superconductors 
2.2.1 Bulk property 

Before extending the above single-band case to a two- 
band one, we briefly review the earlier works by Shul 
et al. 27 ' and Kondo 28 ) for multiband superconductivity. 
The unique property of the multiband is described by 
the following model Hamiltonian: 



1 

T12 



"Hbulk = 



2irn imp N \M L 



/ j / j / j e k^C k ^ a Ck^a 



with 



$2 - Si 
2T r . 



/j.,/j/ — ± kk' 



> 0, A(x) = Si(x)tanhx, B{x) 



5 2 (x)tanhg: erc 



(13) 



4t* 
Si(x) = —ReJ2 

7T * — ' 

n=0 



^(] 



.X 
1 — 

7T 



-n- 2 



-5*2(2;) = ^rlmy^ 

IT 6 *—< 

n=0 



1 

2 

tp(l + n 



.x 

i— 

n 



1 .x 

2~ X n 

Ml 



the first term represents the kinetic energy for the 
±) conduction band and c! (ck^a) is the creation 
(annihilation) operator. The second term represents the 
interaction between electrons with coupling constants 
Vpn' '■ t 1 — I 1 ' f° r the intraband and /i 7^ \j! for the in- 
terband. T c is determined by the following form of the 
linearized gap equation: 



-V++N 0+ 
-V-+N 0+ 



-v+. 
-v__ 

1 



JV _ 
N _ 



A_ 



1 .x 

2 "^ 



A+ 
A_ 



Here, Ti 2 represents the lifetime due to the impurity scat- 
tering. N n is the density of states at the Fermi energy. 
ij)(x) is the digamma function. As a consequence, eq. (10) 
leads to a gap equation for the nonmagnetic scattering 
(T12 is denoted by rf 2 ), 



T 12 lo g7^ = -fA{X) 



f u (x) = f (x) > 0, (14) 



log[(2c^ c )/(7rT c )] l •'• '• (Jii 

Here, A± and A^± are the order parameters and den- 
sity of states at the Fermi energy for the /j, = ± band, 
respectively. 7(~ 0.577) is Euler's constant and we have 
assumed the same cut-off energy uj c (3> T c ) for the two 
bands. T c corresponds to a positive eigenvalue of the 2x2 
matrix on the left-hand side of eq. (17). The multiband 
superconductivity is characterized by whether the inter- 
band interaction is attractive or repulsive. To see this 
point clearly, let us consider a simple case of two iden- 
tical bands, where V++ = V-- < 0, V] = V_+, and 
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N, 



o+ 



as 



Nq- = Nq. The solution of cq. (17) is obtained respectively. By the unitary transformation 

the linearized gap equation is rewritten as 



2e 7 u; c 
T cA = exp 



1 



\V ++ \-V+-)N 



(22) 



T C B = 



2e<uj r 



exp 



1 



(18) 



(\V ++ \ + V+-)N a 
Here, T c a and T cB are the transition temperatures for 

A s = ^±^. A B = ^ T ^, (19) 

respectively When V+- < 0, T cA > T cB (T cB = if 

|1/| | > |V++|). This means that the higher T c is due to 

the attractive intcrband interaction for the order param- 
eters with the same sign (A + A_ > 0). Conversely when 

V-i > 0, the higher T c is obtained for their opposite 

signs (A + A_ < 0). The former and latter are called the 
s ++ -wave and s±-wave states, respectively. 

2.2.2 Impurity effect 

Although electron-phonon and Coulomb interactions 
were taken into account as possible origins of the inter- 
band interaction before, 27 ' 28 ' we here propose inelastic 
(dynamical) impurity scattering as another origin that 
can give rise to T c enhancement in multiband supercon- 
ductors even if it is due to magnetic impurities. For this 
purpose, we extend the work by Fulde et at. to a two- 
band case. 

Then, we reexamine the gap equation in eq. (10) for 
a two-band s-wave superconducting state with order pa- 
rameters, Aju for the /i band. In the absence of impu- 
rities, the superconducting transition temperatures are 
expressed as T c q± for the /j, = ± bands, respectively. 
Here, the interband interaction is not taken into account 

for simplicity (Vi\ = 0). Considering both magnetic and 

nonmagnetic scattering processes, we obtain the follow- 
ing linearized gap equation in the two-band case: 

A + l0g; 7 ' 



8T C 



A_log 



T c o- 



X=0,s 



'12 



Mx)\l + f u {x)>£ 



A^ 



(20) 



where T c is the transition temperature in the presence 
of impurities. The energy difference x between the impu- 
rity singlet ground and multiplet excited states is scaled 
by 2T C [see eq. (13)]. The matrices A^ and A* express 
the impurity intraband and interband scattering contri- 
butions to the self-energies, Ea(iw) and £ w (iw), respec- 
tively: 

it 



A X - 



x x 



(£ = A,w). (21) 

In the following study, it is convenient to use the or- 
der parameters given in eq. (19), where Aa and Ab are 
order parameters for the s ++ -wave and s±-wave states, 



1 , (T c0+ 
2 l0g fe 



8T C] 
= log 

IT 




1 





A(x) 



A A 
A b 



^T c0+ T c0 - I V A b 



(23) 



[2 & \T c0B J 


(i-°0 


+ A(x) 


( A A 
v A B 


ST 
= log 

TT 


f T c 


)U 


^ 


V V ^cOaTcOB 


) 



X=0,s 12 £=A,w 

We can see that T c , which is the superconducting tran- 
sition temperature in the presence of impurities, plays a 
role as an eigenvalue of the matrix on the left-hand side 
in eq. (23). The largest eigenvalue corresponds to the 
highest T c . When T c o+ = T c o- for simplicity, only A(x) 
is left on the left-hand side, and a positive eigenvalue of 
A leads to T c enhancement. 

Next, we consider a case of two identical bands with 
interband interaction (VI) 7^ 0). In this case, the T c val- 
ues arc different for the s ++ -wave and s±-wave states, as 
given by eq. (18). In the absence of impurities, we intro- 
duce transition temperatures T c0 a and T c0B for the s ++ 
and s-|--wave states, respectively. For T c in the presence 
of impurities, the gap equation is expressed as 



(24) 



Here, the contribution of the impurities is expressed by 
A(a;) and it modifies the transition temperatures in the 
bulk. 

In the following part of this paper, we consider only 
two identical bands without the interband scattering 
(V+- = 0) to capture the essence of the T c enhance- 
ment caused by the inelastic scattering impurities, where 
T c oa = T c ob = T c q. We note that it is easy to extend the 
formulation to T c oa 7^ Tcob cases. 

2. 3 Example of magnetic interband scattering for T c en- 
hancement in s±-wave state 
As mentioned above, T c is enhanced by inelastic non- 
magnetic impurity scattering in single-band s-wave su- 
perconductors. In this subsection, we show that magnetic 
interband scattering can also cause T c enhancement in 
multiband cases. For this purpose, we focus on the roles 
of spin-dependent scattering in the s±-wave state, which 
is suggested as one of possible superconducting states 
realized in Fe pnictide superconductors. 12,13 ' As we will 
see, the interband scattering is important for T c enhance- 
ment, which is unique to the multiband and is never seen 
in the single-band case. In fact, neither intraband nor 
elastic impurity scattering can be neglected in real sys- 
tems. These effects cause T c suppression. Whether T c is 
enhanced or not depends on how to balance the pairing 
and depairing effects caused by impurities. 
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Here, we extend the formulation given in §2.1 straight- 
forwardly. The conduction electron part for the p band 
is given by 



^ = E/ dr ^>) e (~ iV )^M r ) 



A u / dr 



^t( r )^i( r ) + VhWVvK^) 



(a* = +,-) 



(25) 



Here, ipfiair) is a field operator of the conduction election 
for the p band. A M is the p band superconducting or- 
der parameter for the s±-wave superconductivity, where 
Ap takes a real value. For spin-dependent intraband and 
intcrband scatterings at the impurities, the interaction 
Hamiltonian is defined by 






R 7 mn uv aa 



X M r 



'V>L(r)Vw(>), 



(26) 



where <r is the Pauli matrix with the three components 
a a (a = 1,2,3). As the corresponding spin exchange, 
M 



mn,uv 



is the scattering matrix element that depends 
on the band components as well. 

We introduce the following 4x4 matrix form of the 
thermal Green's function for the p band: 

G u (t, r, r') = -(T* a (r, r)*t (r', 0)), (27) 

where ^ a {r) and ty* (r) are defined for each band as 



*» 



Vvi( r ) 



*t(r) = ( ^(r) V^(r) ^(r) ^(r) )(28) 

Their Heisenberg representations , JJ^(t*,t) and 'J'Mr,!) 
are given in the same manner as eq. (6). After introducing 
the unperturbed Green's function 

e k p 3 + A M p 2 <7 2 



G^(iwi,fc) 



iui 



^ 



(29) 



as in eq. (7), and combining the two-band forms, we use 
the following 8x8 matrix form of Green's function: 



G (iui,k) 



G+(iwi,k) 

G_{\w u k) 

In particular, for the s±-wave state (A + = — 
it is rewritten as 

e k p 3 + Af 3 /3 2 <7 2 



A_ 



G a (iuJi,k) 



lWj 



(30) 

= A), 

(31) 



M l T" ^ fe -r " 

where f a (a = 1,2,3) is the Pauli matrix for the band 
space. Similarly, the matrix for impurity scattering is 
given by, for instance, 



f> = E E *&»> 



mn a— x,y,z 



U" 



M a 



l T\P 3 (Tl 1 



U v 



M m n Tlfa, 



IP 



MZ mn T\ho 3 . 



(32) 



Here, fi represents the intcrband scattering. For the scat- 
tering matrix U in eq. (32), the self-energy in Fig. 1 is 
given by 

*,. x „ov-v^ 1 1 

E(iwj) 



77- T 

"amp-*- 



2 yy 

„ti iWi - * m iW2 - 5 " 



to 



UZ n G (iuji + iwi - iw 2 , fc)f ° m -(33) 



As in the single-band case, A in eq. (31) changes its 
sign by the U^ n GoU" m transformation. This is a key to 
the T c enhancement by inelastic scattering also in the 
multiband case. 

We apply the above argument to the spin-dependent 
intcrband scattering case such as eq. (32) that enhances 
the T c of the s±-wave superconductivity. In the calcula- 
tion of the self-energy, the scattering matrix U^ n satisfies 

{7«„(Af 3 p2<7 2 )^ ro - HilC„| 2 (Af 3/ ) 2( 7 2 ), 



u a u a = \M a r 

inn nm I mn I ' 



which leads to 

is _ / 1 

Aa - -1 



27rn imp iVo 



Z^ 



A s = 



|Mf 2 | 



1 
1 



(34) 



(35) 



(36) 



in eq. (24). Here, N represents the density of normal 
electron states at the Fermi energy. We have assumed 
here a singlet-singlet configuration for the energy levels 
(m,n — 1,2) of the impurity. Then, we obtain the fol- 
lowing gap equation: 








fo(x) 



A A \ = _^_i___?W A A 

A b ) T c0 g T c0 I A B 



(37) 



where x = (<5 2 — 6i)/(2T c ) and a s represents the strength 
of the spin-dependent impurity scattering defined by 

(38) 



8r c oTf 2 

In the gap equation, fo(x) = —/a {x) + f u {x) > and 
f s (x) = -U(x)-f u (x) > [sec cqs. (14) and (15)]. This 
means that T c is enhanced for the s±-wave state. We 
show the (<5 2 — Si) dependence of T c for various a s values 
in Fig. 2. At Si = J 2 , where the two impurity states are 
degenerate, there is no T c enhancement. When (<5 2 — Si) 
is increased, T c increases accordingly and takes a maxi- 
mum value. In fact, T c depends on two factors competing 
with each other. One is the strength of the attractive in- 
teraction between electrons as derived in the BCS theory. 
The other is the energy region for the attractive interac- 
tion related to the cutoff. The former is intensified by a 
small (<5 2 — Si), while the latter becomes large for a large 
((5 2 — Si). In Fig. 2, one can see the maximum at approx- 
imately (<5 2 — Si) ~ 10T c0 . For (<5 2 — Si) — )■ oo, there is 
no enhancement in T c , since such a higher- lying energy 
level does not contribute to the attractive interaction. 
Thus, the order of (<5 2 — Si) ~ T c is the most appro- 
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Fig. 2. (Color online) (<52—(5i) dependence of T c for various values 
of a s . 



priate for T c enhancement by the interband scattering in 
the s±-wave state. On the other hand, it must be pointed 
out that T c suppression is caused by elastic scattering 
due to the impurity singlet ground state or by intra- 
band scattering [for instance, fi is replaced by unity in 
eq. (32)], which have been neglected here. The most im- 
portant point is the internal structure of impurities that 
intensifies the interband magnetic scattering to overcome 
these pair-breaking effects. 

3. Typical Impurity Interband Scattering for 
T c Enhancement 

In §2.3, we have discussed T c enhancement due to the 
spin-dependent interband scattering in sj--wave states as 
one of the examples, where the details of such impurity 
scattering have been put aside. In this section, we study 
typical impurity scattering for a singlet-multiplet config- 
uration that can give rise to T c enhancement in multi- 
band superconductors. First, we show possible examples 
for the singlet-singlet and next apply the same argument 
to the singlet-multiplet case. In §3.2, we take account of 
not only the interband impurity scattering effect but also 
the intraband impurity scattering effect neglected in the 
previous section. 

3.1 Singlet- singlet configuration 

In addition to the fi type in eq. (32) for the interband 
scattering, there is another type of magnetic scatter- 
ing, M mn T2, due to an orbital moment (spin-independent 
scattering). Since it satisfies 



T2(t3P202)t2 = -(f 3 /3 2 <7 2 ), 



(39) 



T c can be enhanced by the impurity scattering for the 
s± wave. Here, we consider realistic cases for U z nn — 
M^TipzUz and U mn — M mn T2- In practice, we check 
the possible symmetry of Mn for the singlet-singlet con- 
figuration corresponding to each electron scattering type, 
assuming local orbital symmetries of band electrons at 
the impurity site. Although actual bands can include sev- 
eral orbital components, we represent each band by an 
orbital component that mainly contributes to the band 
construction. 

First, let us begin with the spin-dependent case, 



U mn = M mnTl&3 



M z 



(40) 



(V'+t 0+1 V'-t 0-4.) i wnere t denotes transposition 
Uran & given as 

( 1 \ 

0-1 

10 

V -1 / 

To identify the symmetry of electron scattering that de- 
pends on orbital components, we here assume the yz and 
xz orbital types for the /i = + and \x = — bands, respec- 
tively. The wave functions tp y z,tr and ij) X z,a at & n impurity 
are connected to j = 3/2 angular momentum bases ipj z 
(specifically, the Oh Tg point-group bases) by the follow- 
ing unitary transformation as 



rpyzt 

Ipxzf 



1 

V2 



1 

-1 



■03/2 
0-1/2 



(41) 



$vz± \ = J_ I -i -i ) ( 0i/2 

In this new basis set tp j—3/2 = 

003/2 V'l/a ^-1/2 0-3/2)*) it is useful to classify 
the scattering types (dipole, quadrupole, and octupole). 
U^n i n ec l- (40) is transformed as 

/ 



U z 



M z 









i 
















-i 














i 









\ 



/,= 



j—3/2 space 



(42) 

This matrix expression corresponds to the xyz type of 
tensor for j = 3/2. In the local scattering at the impu- 
rity, M^ n has the same symmetry as xyz. In the cubic 
point group, M^ n expresses Oh T 2 octupole coupling (r 3 
for D^h point group). This coupling is realized in the f 2 
configuration (doubly occupied /-electron state). For in- 
stance, it connects the D^h crystal-field ground state 



|.g)= Cl (|4) + |-4))+c 2 |0) (2| Cl | 2 + |c 2 | : 



with the first excited state 
1 



V2 



(|2) 



2)), 



1), 
(43) 



(44) 



in the inelastic impurity scattering U^ n . Here, \M) (M = 
—4 ~ 4) is an eigenstate of J z for the J = 4 angular 
momentum state in the / 2 configuration. Within the two 
crystal-field states, the scattering matrix is expressed by 






M- 
MZ 



12 



\M? 







(45) 



Next, we consider the spin-independent case, U mn = 
M mn T2- In the subspace xp b , we obtain 

/ -i \ 

-i 

i 

\ i / 

Assuming the yz and xz types for the two bands in this 
case as well, U mn is transformed to the j = 3/2 angular 



^ mn — IVl-n 



(46) 



I .. 



= M z nn Tip 3 & 3 . Within the subspace, ip b = 
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momentum basis expression as 



/ 



U„ 



-> Mr, 



-1 














-1 














1 














1 



3=3/2 spa 



(47) 



which expresses the combination of the z dipole and 
z(2z 2 — 3x 2 — 3j/ 2 ) octupole types of j = 3/2 electron scat- 
tering. For the inelastic impurity scattering in this case, 
M mn represents the dipole coupling between two low- 
lying singlet states, e.g., the ground state \g) in eq. (43) 
and the excited state, 

' S(|4)-|-4)), (48) 



V2 y 



given as 



Mn M12 
M 2 i M 22 



= IM12I 



1 



(49) 



Finally, we mention another pair of two-band types, 



(2z 2 



y ) and (x — y ). Applying the above argu- 



ment, we find that the T1P3&3 type of electron scattering 
is reduced to the z(x 2 — y 2 ) octupole type, and the f 2 
type corresponds to the xyz octupole type expressed by 
eq. (42). In Appendix A, we give a different analysis for 

U z . 

mn 

3.2 Singlet- doublet configuration 

Here, we devote ourselves to extending the present 
theory of T c enhancement to typical examples of impu- 
rities with internal degrees of freedom. Orbital degrees 
of freedom of conduction electrons give rise to magnetic 
and nonmagnetic exchange scatterings, both of which 
are considered here in the singlet-multiplet configuration. 
Whether T c is enhanced or suppressed depends on the ra- 
tio of their scattering strengths that determines the signs 
of the two-band superconducting order parameters: they 
are the same (A + A_ > 0) or different (A + A_ < 0). For 
the T c enhancement, the crystal- field ground state must 
be a singlet (a nonmagnetic doublet is also allowed). As 
mentioned for the s±-wave state (A + A_ < 0) in §2.3, 
T c suppression is caused by elastic scattering due to the 
singlet, which is neglected here as well. First, we dis- 
cuss a case of singlet-doublet configuration regarded as 
an S = 1 local pseudo-spin. This S — 1 spin is not a spin 
triplet but a spin and orbitally coupled state, as often 
studied for /-electron systems. The latter can be real- 
ized as an impurity low-lying state in a uniaxial (Dih 
or Deh) crystal field. In the same framework, we tackle 
a more complicated case of singlet-triplet configuration 
discussed in the next subsection. 

For spin and orbitally coupled impurity states, in gen- 
eral, local orbital exchange occurs as well as spin ex- 
change during electron scattering by the impurity mo- 
ment. In the case of an 5* 1 = 1 pseudo-spin for an im- 
purity, a spherical type of exchange interaction is ex- 
pressed by coupling with, for instance, local S c = 3/2 
states, ij} c — (ip3/2 1P1/2 ip-1/2 ^-3/2)*, formed by con- 
duction electrons. This local interaction Hamiltonian is 



given 



by 29 ) 



H, 



i/> ct JsS 1 ■ S c 



rj=l / 



-A 



CF 



3(^) 



(50) 



where a potential (elastic) scattering term is neglected. 
The first and second terms represent dipolar and 
quadrupolar exchanges with the coupling constants, J$ 
and Jq, respectively. The quadrupole operators (Qj_ for 
an impurity; Q't for an electron) are defined as 



{Qn,V= !)••• ,5} 

= \byb z + b z by, b z b x + b x b z , b x by + byb x , 



n2 n2 /rj q 2 

D x D y , (4D Z 



Si - S 2 )/V3}. 



(51) 



The last term in eq. (50) introduces uniaxial (D^ or 
D^h) anisotropy to the impurity states, and Acf (> 0) 
is taken to determine a singlet ground state here. For 
the S l = 1 pseudo-spin, it is sufficient to consider the 
above dipoles and quadrupoles. We do not consider the 
anisotropy of each exchange coupling that usually exists 
in a realistic system, which does not affect the following 
argument. 

On the other hand, we introduce some assumptions to 
electron states as follows. At impurity sites, partial waves 
of conduction electrons are represented by S c = 3/2. This 
can be regarded as the Oh ^s point-group basis in a cubic 
system. In general, actual conduction bands can include 
all four components of S c — 3/2, and their mixing is 
expressed as 



/ 



</> c = V%l> b 



^ 



1p3/2 

^1/2 

lp-1/2 

1p-3/2 



( V'+t \ 



$b 



V'+T 
V'-T 

i'-i. / 

(52) 



where each element in the transformation matrix V is 
given by the overlap of local orbital and band wave func- 
tions such as (S z \[io-) (ji = ±;er =t,|). To examine 
the orbital roles in the superconductivity, we consider 
here the simplest case in which two (±3/2) of the or- 
bital components enter the + band and the other two 
(±1/2) enter the — band. In the above V, we assume that 
(3/2|+ |) = (1/2|- t) = (-1/21- |) = (-3/21+ t) = V 
and that the other matrix elements vanish: 



/ 



V = V Q 



V 






1 














f 














1 


f 












\ 



/ 



(53) 



Such a one to one correspondence clarifies the connec- 
tion between the S c = 3/2 pseudo-spin space and the 
SU(2) spin (£> SU(2) band space, and the scattering ma- 
trices are expressed simply by eqs. (B-9) and (B-I0). Al- 
though this simplification overestimates interband scat- 
tering compared with intraband scattering, it helps us 
examine what types of interband correlations are rele- 
vant to the relative signs of order parameters and how 
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intraband scattering modifies their relative amplitudes. 
For the superconducting order parameters, we con- 
sider both the s++ wave (A+ = A_) and the s± wave 
(A+ = — A_). Since we neglect any correlations between 
the two bands except for the impurity scattering, the 
local correlations directly affect the relative signs of or- 
der parameters together with their amplitudes. This is 
justified when the impurity effect is more relevant than 
any other interband correlation such as an interband 
Coulomb interaction. We respectively express the s ++ - 
wave and s±-wave states as 



A A = Aa/920-2 (s++ wave), 
A B = ABT3P20-2 (s± wave), 



(54) 
(55) 



with the order parameters Aa and Ab , respectively. The 
matrix for the order parameter is given by their combi- 
nation as 



A = A/ 



A T 



(56) 



The relevant impurity scatterings for the singlet- 
doublet configuration are described in Appendix B. First, 
we examine properties of the dipolar and quadrupo- 
lar scatterings separately. Let us start from the dipolar 
(magnetic) scattering case where Js 7^ and Jq = 
in eq. (50). Applying H' in eq. (B-7) to calculate the 
self-energy in Fig. 1, we obtain 

1 1 



£(iwj) = - n imp T 2 ^ V 



^ lWl - S n 1W2 - Sn 

n^n' UJ1LJ2 



V2 



JsVly 

2 I n^ 



oc S~1 r*c I oc fy OC 

Oi_LxqD_ T" 0_LzqO i 



(57) 



where G = G (iw; + iwi — \u 2l k). 5\ (<5 2 ) corresponds to 
the energy level of the impurity singlet ground (doublet 
excited) state (n,n' = 1,2) and 5 2 — Si = A CF . The 
factor (V2) 2 comes from the second-order process S+S- 
or S-S+ between the singlet and doublet states. Then, 
we derive a gap equation. As in eq. (31), the Green's 
function Go consists of three parts: iu>i, ekf>3, and A. The 
£kp3 term disappears after the summation over k. The 
other two parts are transformed by (S+GqS'L+SZGoS'l). 
For A, 

S C + AS C _ + S C _AS C + 

= (5A A - 2A B )/5 2 CT2 + (-2A A - A B )f 3 p 2( T2 

= A A '/0 2 cr2 + AB'r 3 /5 2 <72- ( 58 ) 

This means that the order parameters are transformed 

as 



A'. 



A 



v\ 



A 



5 -2 
-2 -1 



(59) 
Similarly, the iw; part is 

SfywiSl + SliuiSl = ioji{5 - 2f 3 ). (60) 

Since f 3 is the Pauli matrix for the ± band space, eq. (21) 




k = 0.25 



k = 0.75 
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Fig. 3. (Color online) x [= Acf/(2T c )] dependence of the highest 
eigenvalue of the matrix Ti2A(x) for T c . The plot is shown for 
fixed k = Jr>/Jg values. 



is obtained as 



X s , = 



for X = S and £ = w. Then, eq. (22) yields 



(61) 



(62) 



Both A^ and A^ can be divided into interband and in- 
traband scattering parts as follows: 



A A — 



A 







3 



3 
3 



2 -2 
-2 2 



(63) 



If the intraband contribution, the second matrix on the 
right-hand side of each equation, is removed, the inter- 
band contribution results in the T c enhancement for the 
s± state, as discussed in §2.3. 

In the same manner, the above argument is applied 
to the quadrupolar (nonmagnetic) scattering case where 
Js = and Jq 7^ in eq. (50). For the A and iw; parts 
in the self-energy, 



Q C + AQ C L + Q1AQI = -\2A K p 2 a 2 + 12A B T 3 p 2 <J2, 



Q c + iuiQ c _ + Q c _icJiQ c + = 12ioj h 
respectively, and we obtain 

19. = 



a2 = 



-12 
12 



12 
12 



(64) 
(65) 

(66) 



These lead to the T c enhancement for the s ++ -wave state. 
Next, we see the more generic treatment combining 
both magnetic (S) and nonmagnetic (Q) scattering terms 
as 

fiaA(a:) - cosC /a(z)A| + /„(a?)Af 

+ sinC[/A(ar)Ag + / w (x)A« 
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Fig. 4. Ratio /o//s as a function of x = Acf/(2T c ). 



/ [-5/ s (x)cosC 

+ 12/o (x) sin C] 

\ 2/ s (x)cos< 



2/ s (x)cosC \ 

/ 



{[3/ (x)-2/ s (x)]cosC 
-12/ s (x)sinC} 



in eq. (24), where 



T12 = 



T s T Q 

'12 T 12 



J 2 

tan£ = —j = k, 



(67) 



(68) 



and the lifetime rf 2 {j\2J is introduced for the magnetic 
(nonmagnetic) scattering. The highest eigenvalue of f 12 A 
determines the quantity of 



8T C _ T c 

Tl2l0g^-, 

71" JcO 



(69) 



and its x dependence is shown in Fig. 3. The T c en- 
hancement is obtained for its positive value, which holds 
at x > 6.2 for k = and at x > 5.8 for k = 1.0. In the 
presence of only nonmagnetic scattering for k = 00, T c is 
enhanced in the entire x > range. The T c suppression, 
found for a small x or k ~ 0.25, is due to the compe- 
tition between the magnetic and nonmagnetic exchange 
scatterings and to the depairing effect by the intraband 
scattering. Since /o// s < 1, the condition for a positive 
value of eq. (69) is given by 

(-5/. + 12k/ )(3/o - 2/ s - 12«/ B ) - (2/ s ) 2 < 0, (70) 

the solution of which determines the minimum /o// s as 
a function of k: 



1 
24k" 



[(48k 2 + 8k + 5) 



-v/(48k 2 + 8k + 5) 2 - 96k(10k + 1) . (71) 



The x dependence of /o// s and the k dependence of 
(/o//s)min are shown in Figs. 4 and 5, respectively. At 
both k = 0.03 and k = 0.75, for instance, eq. (71) takes 
almost the same value, ~ 0.51, which gives the minimum 
x ~ 14.5 for T c enhancement, as shown in Fig. 3. The 
maximum at k = 0.25 in Fig. 5 indicates that there is no 
T c enhancement for any finite crystal-field level splitting, 
implying that T c is always suppressed by the competi- 



Fig. 5. k (= Jn/Jg) dependence of the minimum of the ratio 
fo/fs for Tc enhancement. 



tion between the magnetic and nonmagnetic scattering 
effects. We also find that k determines the relative signs 
of the two order parameters: A + A_ < for k < 0.25 
and A + A_ > for k > 0.25. The combination of Aa 
(s++ wave) and Ab (s± wave) is caused by the intraband 
scattering effect in the — band, due to the off-diagonal 
elements in A(x), which leads to |A_| < |A + | in the 
vicinity of T c . At k = 0.25, A_ = (A A = A B ) means 
that only one band (+ band) is superconducting with 
T c suppression. 

The above argument is based on a rather artificial 
assumption about the local band character in eq. (53). 
More generic treatment of V in eq. (52) provides us with 
various scattering effects on T c . This point is considered 
for the singlet-triplet configuration discussed below. 

3.3 Singlet-triplet configuration 

The singlet-triplet configuration is realized for a non- 
Kramers ion in an Oh crystal field environment like the 
p r 3+ or tj4+ j2 i ow _iy m g states in heavy-fermion materi- 
als. For the strong spin-orbit coupling, the most relevant 
local /-electron states are described by the j = 5/2 an- 
gular momentum. Then, we consider only the exchange 
coupling between the impurity states and the j = 5/2 
electrons hybridized with conduction bands. We assume 
here that the Tg and T^ partial waves (see Appendix C) 
arc transferred independently to the + and — bands, 
respectively. In terms of eq. (C-6) for both r§ and V-?, 
tPm = (V'mit VVni ipm 2 t "0m 2 i ipmrf Vv^l)* is com- 
bined with ip b = (V>+t ip+i ip-f "0-i)* f° r the bands. 
This is expressed by ip m — Vipb'- 
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V 
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U+Mlt 


V+Vi^ 
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v+vu 

v + u 2 -\ 


V + V2f 











W+«2i 




V+U21 






V- 















V- 



\ 



(72) 



/ 



Here, v± represents the hybridization amplitude of the 
/-orbitals and the jjl = ± band, respectively, at the im- 
purity sites; 

Wit = ( m * t l+t), ^t = ( m *tl+l), 

Vii = {m i i\+ J [), un = (mil\+l) (73) 
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represent the impurity site overlaps of the wave func- 
tions iprni.a- (« = 1,2) and ip+,o-- We only simplify the 
connection between the 771,3 orbital and the — band. 

In Appendix C, the relevant impurity scatterings are 
described for the singlet-triplet configuration Ti©r (T = 
r4,rg). For W in eq. (C-23), the self-energy in Fig. 1 is 
obtained as 

1 1 



E (kji) = - 77 imp T 



E 



E 

X=S,Q 



2 E 

n^n' UJ1UJ2 

.2 1 y^ 
k 



(J T xY 



iwi - 5 n 1U2 - 5 n > 

{vc.y s^i vc.r 
A z ' U A Z 



rc,r / 



^C.l si v c, 

A^_ Lj[)A_ 



c.r 



7c,rA x^c,r 



[G = G (iuji +iuj 1 -iw 2 ,fc)]. 



(74) 



Here, <5i (62) corresponds to the energy level of the impu- 
rity singlet ground (triplet excited) state (77,77' = 1,2). 
For the gap equation in eq. (24), the transformation to 
A in eq. (56), 

X c z - r AX c z ' r + -(X c /AX c _: r + l!' r Alf), (75) 



leads to 

3 2 
Al = (-l)A^ = g E( Ul t u a ~ v tt v ti ) 

i—1 



-1 
1 



(76) 



in both r = T4 and r = T^ cases. For the derivation, we 
use UifUii = u**.u*< and Uj-j-Uij, = v*-t-V*\ since the m^ and 
m^ local electrons are the time reversal partners. In the 
same manner, for both T4 and T5, 

Xl> T \uiX c / + ^(XfiunX^ + X c J r iuJiXf) (77) 



gives 



»=i <T=t4 



+ V* a V ia ) 



1 

1 



(78) 



Equation (76) implies the competition between the mag- 
netic (S) and nonmagnetic (Q) scattering effects for 
T c enhancement. Which has the higher T c , the s++ 
wave or the s± wave, depends on whether the sign of 
(uifUii — VifVii) is positive or negative, respectively. In 
realistic systems, both magnetic and nonmagnetic scat- 
terings coexist. In the present case, we can calculate the 
coupling constants, Jg and Jq, based on the Anderson 
model including the j — 5/2 electron exchange scattering 
due to a single impurity in the f 2 configuration, 30 ' and 
obtain the ratio as | Jq\/\ Jg\ — 1/3. 31 ) This means that 
the magnetic interband scattering dominates the T c en- 
hancement for the singlct-triplct configuration. A more 
detailed discussion is given below. 

4. Pr Impurity Effect in LaOs4Sbi2 Supercon- 
ductor 

The Pr 3+ f 2 configuration is a good candidate for 
raising T c if Pr can be embedded in a multiband su- 



perconductor. In fact, for the skutterudite superconduc- 
tor Lai_ a; Pr 2: Os4Sbi2, the Pr singlet-triplet configura- 
tion may be relevant to T c enhancement in LaC^Sb]^. 25 ' 
Here, we show an attempt to apply the above argument 
in this case. The most intriguing feature of the skut- 
terudite compounds is that each rare-earth ion is located 
at the center of the pnictogen cage (Sbi.2) having the 
a u (xyz) and t u (x,y,z) molecular orbitals. It is consid- 
ered that the Pr /-electron states hybridize with the 
conduction bands via these orbitals. For a strong spin- 
orbit coupling, the a u electrons have the Oh ^7 symmetry 
and transfer directly to the 7773 electron state named in 
eq. (C-6). On the other hand, the t u electrons with the 
Oh, Fg symmetry mix with both mi and 7772 states as 32 ' 



f 



|mi,t}++ -j={\x,i)-\\y,i)), 
|mi,i) ++-— (| x ,t)+%,t», 



|m 2 ,t) *+ 



|ra 2 ,T) <-► 



V3l 



V3 



\/2>,t) 

V2\z,l) 



V2 



V2 



(\x,i)+i\y,i)) 



(M>-%,t» 



(79) 



(80) 



, (81) 



(82) 



We here consider two conduction bands, one of which is 
a u -dominant and the other is i u -dominant. We assume 
that both are combined with each other only through 
interband electron scattering. In the present case, the 
most relevant is the Pr impurity scattering due to the 
hybridization effect with the bands. Then, we can use 
the transformation in eq. (72) for mixing /-electron 
states with the t„-dominant band (+ band) and the a u - 
dominant band (— band). Choosing the local t u (x,y,z) 
component that mainly contributes to the + band and 
taking their onsite overlaps arbitrarily as 

(x\+) : (y\+) : (z\+) = sin 9 cos <fi '■ sin 9 sin (f> : cos9, 

(83) 



we have 



(1/V2)sin0 e" 



wit = U U = °i 

7; lt = (l/V2)sin0e^, v u 

7/,2-f = u 2 _i = (v2/3) cos 9, 

v n = -(1/V6)sm9 e _i ^, v 2i = (1/V6)sm9 c i,#> (84) 
Substituting them in eqs. (76) and (78), we obtain 

Ai^-DA^if- 1 ;' 



Af = Ag = - 



1 
1 



(85) 
(86) 



In the gap equation (24), we have 
1 



A(x) 



A T S 
* T 12 



4r 



12 



fo(x) 

-f s (x) 

-/.(*) 

fo(x) 



(87) 
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In the real system, we must consider the Th symme- 
try that combines the Pr Oh T4 and T$ triplet states 



as 



31,33,34) 



|ri 2) r ? ) = v / i-rf 2 |r 5 r/> + d|r 4 7 ? > (77 = +,0,-), (88) 

where d (0 < d < \j\pi) represents the deviation from 
the Oh symmetry In the above calculation of the self- 
energy, we usually have additional terms, namely, 



X c > l4 AX', 



Xb 



2 y 



c,r< 



AX 



c,r 5 



X C ' F4 AX_ 



c,r 5 N 



+ (r 4 ^r 5 ) (x c = s , q ), 



(89) 



which vanish in this case. They also vanish when A is 
replaced with ilo. If the a u band hybridizes with the /- 
orbitals more strongly than the t u band, which means 
u_ 3> w+ here, the impurity interband scattering en- 
hances T c for a finite crystal-field splitting x [= (S 2 — 
S 1 )/(2T C )}, satisfying either f /f s > t? 2 /t? 2 or / // s > 
Ti 2 /ti2- For the former, the T c enhancement is possible 
in the s++-wave state, and a higher T c can be realized as 
the magnetic scattering becomes more dominant. In §3.3, 
we have mentioned that | Jq|/| Js| = 1/3 (t^/tH = 1/9) 
is satisfied for Oh- It also holds for 2\, 31 ) so that the 
gap equation chooses the s++ wave for the T c enhance- 
ment. One may think that the quadrupolar scattering is 
the most relevant if d is small in eq. (88), since the Fi- 
rs interchange dominates the exchange scattering. For 
the T§ scattering type, however, the octupolar scattering 
cannot be neglected in the multiorbital exchange owing 
to the hybridization of /-electrons with the conduction 
bands. Thus, the magnetic scattering can contribute to 
the T c enhancement for the Ti-r^ configuration in a 
multiband system. Our result indicates that the s ++ - 
wave state is favorable for the Lai_ :E Pr a ;Os4Sbi2 super- 
conductivity if the multiband picture is applicable and 
the Pr or Sbi2 site local orbital symmetries reflect in the 
bands. 

Finally, we mention the effect of intraband scattering 
neglected in the above argument where v_ ^$> v + is as- 
sumed for the hybridization amplitudes. If the m\-m 2 
scattering terms in eqs. (CT7) and (CT8) are consid- 
ered in calculating the self-energy for the T5 nonmagnetic 
scattering type, a correction term, 

„,2 



C v —fo(x) ( 1 ■ 



'12 



^v Oil 



(90) 



is added on the right-hand side in eq. (87). It assists the 
T c enhancement, which resembles the effect of inelastic 
nonmagnetic scattering impurities in single-band s-wave 
superconductors. 

5. Conclusion 

We have studied inelastic (dynamical) impurity scat- 
tering effects on T c enhancement in the two-band su- 
perconducting states, s++ wave and s± wave. The key 
is to check the sign change of the order parameters in 
the self-energy corresponding to the second order of the 
scattering process. We solve a gap equation, where both 



the s ++ wave and s± wave are combined by impurity 
scattering, and find out the possible atomic structure of 
the impurity that increases T c . For the T c enhancement, 
the crystal-field ground state must be a singlet or a non- 
magnetic doublet. If not, magnetic impurities will always 
cause T c suppression. For the s± wave, it is necessary for 
elastic (nonmagnetic) scattering to be relatively small. In 
the singlet-singlet configuration, we can easily determine 
what type of magnetic interband scattering contributes 
to pairing interaction. This simple analysis is very useful 
for determining a scattering type for the T c enhancement 
among multiorbital interaction terms when we consider 
such a complicated atomic structure as the /" configu- 
ration. Here, we have discussed the singlet-doublet and 
singlet-triplet configurations. We find that the T c of ei- 
ther the s++-wave or s±-wave state can be enhanced by 
the dynamical magnetic or nonmagnetic impurity inter- 
band scattering for a larger crystal-field splitting of im- 
purity ground and excited states, while T c is suppressed if 
both scattering strengths are comparable. Whether T c is 
enhanced or not also depends on how the local electron 
states hybridize with the two bands. We show a case of 
T c enhancement in the s++ wave by magnetic interband 
scattering due to the singlet-triplet configuration. This 
result may give useful information on a multiband pic- 
ture of the Lai_ 2; Pr a ;Os4Sbi2 superconductivity. In fact, 
the connection between the La-rich and Pr-rich super- 
conductors is left to be clarified as well as the symmetry 
of the latter order parameter. 35 ' 

In the above argument, we have assumed weak scat- 
tering impurities. Strictly speaking, dynamical scatter- 
ing effects should be investigated as a Kondo problem. 
It was shown theoretically in normal metallic cases that 
the impurity exchange scattering strength is renormal- 
ized to be weaker with the decrease in temperature since 
the crystal-field singlet competes with the Kondo-singlet 
formation. 29, 36 ~ 38 ) This holds even for a small crystal- 
field splitting of the singlet ground and excited multiplct 
states. Thus, our results shown here are valid for the 
T c enhancement due to the impurity scattering. 

Throughout the paper, we focus on the hybridization 
between local orbitals and conduction electron states. For 
most of the /-electron systems, this may be more rele- 
vant than the admixture of orbitals formed by intersite 
electron hopping that we have neglected here. Our treat- 
ment is more practical since the latter contribution can 
be included effectively in the hybridization. 

The possibility of the s±-wave state has been the sub- 
ject of debate for high-T c superconductors with FeAs lay- 
ers. It is pointed out that the Fe d-orbitals contribute 
to the disconnected Fermi surfaces and the orbital de- 
grees of freedom play an important role in pairing in- 
teraction. 8 ' The interband Coulomb interactions may be 
relevant for the high T c of this multiband superconduc- 
tivity, 28 ' 39 > 40 ) the roles of which are analogous to those of 
the local correlations in the /-electron systems we have 
considered here. 

Acknowledgement 

This work is supported by a Grant-in-Aid for Scien- 
tific Research (No. 20540353) from the Japan Society for 



12 J. Phys. Soc. Jpn. 



Full Paper 



the Promotion of Science. One of the authors (H. K.) is 
supported by a Grant-in- Aid for Scientific Research on 
Innovative Areas "Heavy Electrons" (No. 20102008) from 
The Ministry of Education, Culture, Sports, Science and 
Technology, Japan. 



Appendix A: Connection between the s-|--Wave 
and Conventional s-Wave Super- 
conductors 



We attempt to apply the unitary transformation 
M -1 1 









rise to the T c enhancement if the electron scattering in- 
terchanges the impurity low-lying states with a finite 
crystal-field splitting. 

Appendix B: Exchange Matrices for Singlet- 
Doublet Configuration 

At low temperatures, the most relevant terms in 
cq. (50) are electron scatterings with interchange be- 
tween the local ground and excited states given by 

Js [SlS c x + SlS c y ] + J Q [Q\Q\ + Q\Q C 2 ] . (B-1) 

(A-l) Using S± = S x ±iS y and Q± — ±iQi + Q2, it is rewritten 



to the band basis, where M is a natural number. Using 
this new basis set, ip M = (ip M f V'mj. tp-Mf tp-Mif, 
U^ nn in eq. (40) is transformed as 

/ 1 
0-10 



l) z 



-M z 





V o 



-1 






-1 








1 ) 



M space 

(A-2) 



If the number M is regarded as an orbital component, 
this matrix indicates a quadrupole type of scattering. In 
fact, ij>M=i corresponds directly to tpj =3 , 2 in §3.1. In 
this case, the matrix in eq. (A-2) is expressed by (2jf — 

Jx Jy)' 

On the other hand, for the conventional single-band s- 
wave superconductivity, the pairing Hamiltonian is given 

by 



h a = - E A ( c U c -h + c-HCk-t 



(A-3) 



After applying spherical expansion to the operators as 

/g^-jl/2 
Cka =X/ kR Y LM{tt k )c kLM a, (A-4) 

LM 

where R is the radius of the system, we rewrite the 
Hamiltonian (A-3) as 

h a =-y; e m) Ma 

k L,M>0 
X ( c feLMt c feL,-Aa + c kL.,-MiC kLM \j ■ (A-5) 

This orbital component M can correspond to M in 
eq. (A-l) directly. The pairing interaction works between 
the electrons with the different orbitals denoted by ±M. 
Let us fix L and restrict the orbitals to a pair of ±M. If 
ipMa is transformed to ip fl(7 using eq. (A-l), eq. (A-5) is 
reduced to 

^a = - E E(-!) A VA (4„t c L + c */4<Vt) ■ 



V=± k 



(A-6) 



It represents s±-wave pairing in the two orbitals con- 
nected to the \x bands. 

Thus the spin-dependent intcrband scattering in 
eq. (40) in a two-band s-t-wave superconductor can be 
mapped to the quadrupole type of scattering in eq. (A-2) 
in a single-band s-wave superconductor. Both cases give 



^-[SlSl + SlS^+^-iQlQl. 



(B-2) 



By denoting S l z = (the impurity ground state) and 
S l z = 1,-1 (the impurity excited states) by |1), |2), and 
|3), respectively, the (S±) mn and {Q l ±) mn (m, n = 1, 2, 3) 
matrix expressions are given as 



si 



Ql 





( ° 





V2 \ 


s [ _! = 


V2 










V o 





o J 




( ° 





-V2 


Q v l = 


V2 










\ o 









(B-3) 



(B-4) 



For the S c = 3/2 states, the operators S'L = S_ and 
Q c + = Q_ are expressed by 



/ ^/3 



^S c + ip c = i/> 



c - „/, c t 






Vo 










\ 



v> c 



J 



/ 2\/3 



^Qlip c = V 



cf 















V o 









-2V3 





(B-5) 



t/> c . (B-6) 



Following eq. (26), we examine the low-temperature 
physics using the impurity interaction Hamiltonian 



W = E J2J2 dra\ m a in S{r - it, 



R~ mn 



Jx 



x *t (r) _* (xl imn X c _ + X\ mn Xl) *(r). 

(B-7) 

Here, ^ is an eight-dimensional vector for conduction 
electrons that is obtained by extending \J3 b in eq. (52) to 
the particle-hole space as 



* 






*T = ( * 



*I 



(B- 



where \I/ M (/z = +,— ) is defined in cq. (28). Through 
V in eq. (53) in this extended space (Pauli matrices a 
for spin, p for particle hole, and f for band), the above 
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electron scattering matrices are rewritten as 



S+ = — ti(/3 3 <ji -io- 2 ) 



■^i 1 -h)(P3&i +102), 



q; = V3(if 2 )( 



0"! - 1P302) 



(B-9) 
(B-10) 

where \Vo\ has been omitted. One can see that S± con- 
sists of the interband scattering n and intraband scatter- 
ing (1— T3) terms. The former is comparable to (U x — \U y ) 
in eq. (32). 

Appendix C: Exchange Matrices for Singlet- 
Triplet Configuration 

On the basis of the J = 4 total angular momentum, 
the Oh singlet and triplet states are given by 34 ) 



|ri> 



(|4> 



^(14) 



13) 



4)) 
•3>- 



12 

|r 4 +) = 
|r 4 o) = 

|r 4 -> = 

|r 5 +> = 
|r 5 o> = 



|rr,-) = - A /g|-3) 



4», 



(C-l) 



(C-2) 



(|2)- 



-1 



2», 



(C-3) 



> 



For the j = 5/2 local electrons, they are classified into 
the Oh symmetric states as 32 ) 



^8,3/2) - 

l r 8,l/2) = 
l r 8,-l/2> 



- 8,-3/2/ 



■7,1/2/ 



7,-1/2/ 



6 
|l/2), 



13/2) 



1/2), 



^1-3/2) 



5/2), 



>/2), 



(C-4) 



||-3/2)-^|5/2>, 

1, 



(C-5) 



|3/2)- A /-|-5/2) 



In the following argument, we change their notations as 



\-3/2 -> ( m i,t), r 83 / 2 



8,1/2 



(m 2 ,t), 
(m 3 ,t), 



8,-1/2 



(mi, 4.), 

(m 2 ,4-); 
(m 3 ,4-). 



(C-6) 



17,1/2 -* \jn,3, 1 j, 17,-1/2 

In the same manner as that in eq. (B-3) for an impurity, 
the magnetic coupling between the Ti singlet ground (de- 
noted by |1)) and T (= r 4 ,r 5 ) triplet (+, and — states 
denoted by |2), |3) and |4), respectively) excited states 
is expressed by the following matrices: 31 -' 

/ 1 ^ 



10 

\ / 



ij ~ — 



(C-7) 



ij I — ij _ 



(C-8) 



Ql = 



/ v/2 \ 

ii = -V2 

0' 

V 000/ 

For nonmagnetic coupling, 

/ -i \ 



i 

\ / 

/ -y/2 \ 
-V2 


V 00 / 

We note that the relative signs of the off-diagonal matrix 
elements are different between S± (Q±) for the pseudo- 
quartet and S± (Q±) for the S l = 1 pseudo-spin given 
in eq. (B-3) [in eq. (B-4)]. The singlet-triplet interchange 
processes are mediated by the corresponding Tg and T7 
electron scatterings as follows: 31 -* 



Ql = Q- 



(C-9) 



(C-10) 



■fptn&tPm = 2 W'mit^t + V^Vw) + Hc -> ( C ' n ) 

+ ^(V-U^3t-^ 3 ^ ro2 t), (C-12) 

V'L^V'm = ^(V^^msT + V^+Vw) + B.C., (C-13) 

+ ^(-^ mi ^m 3 t + ^ m3 ^ m ^), (C-14) 



ll>\n<Fz i, $m = -^(V'mjt^t + V^Vw) + B.C., 



(C-15) 



^L<F+^ m = ^rm^msi + i>i l3t ^ mi i) 



2 V ™m 2 |''" l 3T V m3 ,p 



+ ^(-V-L, V-mst - *l>Ll*l>mrt), (C-16) 



^InQ^^m = ~i j 2 ^m 2 tV , m 3 t + V^Vw) 

+2\/5(^ lt V' m2 T " V'L+VW)} + H - c - 

(C-17) 

+ — (^m 1 4.V , m 3 t + ^ma^it) 

In §3.3, we assume v_ 3> w+ for the hybridization with 
the conduction band in eq. (72), i.e., we neglect the m\- 
m 2 orbital scattering terms in eqs. (C-17) and (C-18), 
whose contribution is small in the band scattering. Using 
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the transformation in cq. (72), the above multiorbital 
scattering operators are reduced to interband scattering 
operators. For the T4 type, 



i>\^b = i>\\ 



( 






wit 
V v n 






wit 



ris T + *iP b = iPt- 
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"it 

"it 









U14. 
-\/3v2t 



7 U 



\ 






>/3«24. 







</>(, 



(C-19) 



2. 



'it 
'it 



-v 


(C-20) 



</>6 



/ 



^^ 



^ 



/ 



-wit 
V -vn 






-wit 

-W14. 



/ 



^#^=1^ 







i>H U14. 

\ — v3«2t — v / 3i'2t 



"it 

"it 






-y/3v^, 

■y/3v£± 












(C-21) 



u it 

fit 




(C-22 



^6 



where (v+t>_), which appears as the common factor in 
each matrix element, has been omitted. For the Tc, type, 
s^ 5 is obtained by switching the indices as 1 -B- 2 in the 
s^ 4 matrix, and s + s by replacing V3 — > — v3 in addi- 
tion to switching 1 -R- 2 in the s + 4 matrix; q^ 5 and q + 5 
are obtained from q^ 4 and g + 4 , respectively, in the same 
manner. We note that s^ 4 in eq. (C-19) is comparable to 
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12 
13 

15 

16 

17 
18 

19 



C 7 ™™ in eq. (40) if u lt 
,r 4 



-W14, 



'it 



are taken; q], 4 in eq. (C-21) is comparable to U mn 



'it 



'"U 



* 1 and v\ a = 
bl< 
and via- 



111 
0. These 



cq. (46) for wit = U U 

scatterings, whichever is magnetic or nonmagnetic, con- 
tribute to T c enhancement in the s±-wave state. As dis- 
cussed in §3.3, whether T c is enhanced or suppressed by 
the interband scattering for the s± wave and also for the 
s++ wave depends on the details of electron transfer, Ui a 
and Via, between the /-orbitals and conduction bands at 
an impurity. As in cq. (B-7), we introduce the impurity 
interaction Hamiltonian for the singlet-triplet configura- 
tion for the r (= r^Ts) scattering, 



X=S,Q R-, mn 



dra\ m a jn S{r 



Ry) 



&{r)J r 
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X X 

I z .mn 



c.I' 



>' 
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X X 
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where the band electron scattering operators 



jc,r 



Q 



c,r 



Q 



*(r), 

(C-23) 

(C-24) 



are rewritten on the basis of the eight-dimensional vector 
* in eq. (B-8). 
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